Generation of quantum discord between ionic qubits via noisy processes 
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Quantum systems in mixed states can be unentangled and yet still correlated in a way that is not possible for 
classical systems. These correlations can be quantified by the quantum discord and might provide a resource 
for certain mixed-state quantum information processing tasks. Here we report on the generation of discordant 
states of two trapped atomic ions via Markovian decoherence processes. While entanglement is strictly non- 
increasing under such operations, discord can be generated in various forms. Firstly we show that, starting from 
two classically correlated qubits, it is possible to generate discord by applying decoherence to just one of them. 
Secondly, even when starting with completely uncorrelated systems, we show that discord can be generated via 
classically correlated decoherence processes. Finally, the Werner states are created. The generated states can be 
used as a resource state for quantum information transmission and could be readily extended to more ions. 



Over the last decade improvements in experimental con- 
trol over quantum systems have enabled scientists to gener- 
ate non-classical states of light and matter. Entanglement be- 
tween quantum systems [ 1 1 has been at the forefront of this 
research largely due to the key role it plays in Bell inequal- 
ities [2 1 and as a resource for certain information processing 
tasks @. However, entanglement is not the only kind of non- 
classical correlation that can exist between systems. Discord 
was introduced as a measure of these other correlations in bi- 
partite systems, which becomes distinguishable from entan- 
glement for mixed states [4-7 1. There are mixed states with 
discord yet no entanglement and these are the focus of the 
present work. 

There is growing evidence that, like entanglement, discord 
can be viewed as a resource for information processing tasks 
ll6ll8 UT2l . Recently, quantum states of light have been gen- 
erated which belong to the regime where there is discord but 
no entanglement ||8l [TUl [13l [14ll . Following extension of the 
discord concept to continuous variable systems |[T5l[T6l . ex- 
periments have also investigated Gaussian discord dynamics 
under various decoherence channels [14|. Here we are con- 
cerned with discrete variable systems (qubits). 

This paper is organised as follows: After briefly reviewing 
the definition of discord we demonstrate how operations on 
one qubit can generate discord between two qubits, via a sim- 
ple example. Secondly, we review the correlation rank as a 
way to assess the nature of correlations in quantum states and 
then show how discordant states with any rank can be gener- 
ated via noisy processes. Finally, we present the first genera- 
tion of the Werner states between atomic qubits. 

Two quantum systems A and B have discord, when con- 
sidering measurements on system A, if and only if their state 
cannot be written in the form pab- 2/ Pil^f )('Af I ®pf> where 
(i//i\i/fj)=dij, pf are density matrices of qubit B and p, are prob- 
abilities. For a state pab a von Neumann measurement Q of 
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FIG. 1: Potential energy level scheme of the 40 Ca + ion showing 
the relevant Zeeman sublevels. A qubit is encoded in the states 
S \ji(nij= + 1/2)=|0) and £>5/2(m ; = + 3/2)=|l), where m, is the mag- 
netic quantum number, ancillary levels \d) and \d) are temporarily 
employed for state preparation, as described in the text. The qubit 
is at optical frequencies (729 nm), the re-pump is infrared (854 nm) 
and the decay channel is UV (397 nm). Adjacent Zeeman states in 
the S and D manifolds are shifted by about 10 MHz. 



A with eigenvectors n,=|t/r,)(i/r,| will leave the total state un- 
changed, i.e. Yji^iPAB^-PAB- The discord D of a bipartite 
system is quantified as the difference between two definitions 
of the mutual information, / and J, i.e. D=I—J. The first def- 
inition (/) captures the total correlations of the density matrix 
Pab by the difference in entropy of systems when taken in- 
dividually and when taken together I(pab)-S (pa) + S(pb) - 
S(pab), where S is the von Neumann entropy |3| and Pa(.Pb) 
is the reduced density matrix of system A (B). The second 
definition (J) captures the classical correlations and can be in- 
terpreted as information gain about one subsystem as a result 
of a measurement on the other. J(pab)=S (pb)S (pb\a), where 
S (pb\a) 1S the entropy of B after a measurement of A (with 



unknown result) and J is maximised over all von Neumann 
measurements of A, see Q. Discord can be asymmetric with 
respect to exchange of the two systems since it quantifies the 
extent to which measurements on one system affects the total 
system for an independent observer. For the discord that con- 
siders measurements on system A (B) we will use the label 
D A (D B ). 

Consider the separable two-qubit mixed state 
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(I+>a<+I® !+>*<+! + I->a<-I«I->*<-D, (1) 



where |±)=(|0)+|1))/ V2, Although manifestly correlated, 
state piis fully classical since it is diagonal in the local orthog- 
onal |±) ® |±) basis. The discord is correspondingly zero in ei- 
ther direction: a von Neumann measurement {II ± = |±)(+|) 
of either qubit leaves the density matrix unchanged. Sur- 
prisingly, discord can be generated by applying an amplitude 
damping process that acts only on one of the qubits in state pi 
lfT71fT8l. 

Amplitude damping of a single qubit can be described by 
the quantum map e' if/ (p)=£ , op£' ( ! ) +£ , ip£ , j with Kraus operators 

£o=|0><0|+ VT^pi'lXll and £,= VpI0><1|. Forp=0 (p=l) the 
qubit undergoes zero (complete) amplitude damping. Con- 
sider now the case where after preparation of pi qubit B un- 
dergoes this process, i.e. it interacts with a dissipative Marko- 
vian bath that causes the excited state |1) to decay to the 
ground state |0). Throughout the damping process the state 
is of the form e ad (pi)=\(\+) A (+\ ® t +b + |-)a<-|®t_ b ), where 
T ± B=e' ad {\±)(±\) are generic density matrices representing the 
state of qubit B. The effect of the damping process is to re- 
duce the distinguishability of t±b, which are initially orthog- 
onal. Subsequently, qubit A becomes correlated with non- 
orthogonal states of qubit B. For 0<p<l there is no von Neu- 
mann measurement of B after which the original state, e a d(pi), 
is recovered and therefore the discord Db>0- 

Experiments are carried out using a system of two 40 Ca + 
ions in a linear Paul trap. A qubit is encoded in an S 1/2 
ground and a D5/2 metastable state, which are coupled by a 
narrow quadrupole transition (see Fig. IT). Correlations are 
established between the state of the valence electrons in two 
ions. To generate the complex correlated mixed states addi- 
tional electronic levels are employed to both temporarily store 
information and remove entropy from the system. 

Firstly, the initial classically-correlated state p\ (Eq. 1) is 
prepared. Using a standard optical pumping technique, both 
qubits are initialised into |0) (Fig. [T). The two qubit unitary 
operation MS (6-n / 4)= exp(-/#<xjcrj) is applied, where o~ x is 
a Pauli spin- 1/2 operator, generating a maximally entangled 
state (|00)-/|11))/ V2 |19|. This operation is realised using a 
narrow-band bichromatic laser pulse with frequency compo- 
nents close to the first red and blue axial centre-of-mass vibra- 
tional sidebands of the two-ion crystal, as described in detail 
in [20 1 . A 20 kHz symmetric detuning from the sidebands 
is employed, entangling the qubits in 50 /-is. Next, a series 
of resonant 729 nm pulses moves this superposition into the 
state {\dd)-i\aa))l y/2. An 854 nm laser pulse then couples 
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FIG. 2: Amplitude damping results e a dipi) derived from experimen- 
tally reconstructed density matrices, (a) Discord D A (ideal dashed 
blue line, data squares) and D B (ideal green line, data triangles) as a 
function of amplitude damping probability p (see text), (b) Ideal 
trajectories T +B =e' ttd (\+X+D (blue solid line) and r_ s =e^(|->(-|) 
(red dashed line) on the X-Z plane of the Bloch sphere. Open 
shapes with error bars show experimental results derived from the 
reduced state of qubit B in Eq. (1) after projecting qubit A, in 
the reconstructed density matrix, into either |+> (triangles) or |— ) 
(squares) i.e. Tr / 4(|±>(±|e flrf (p 1 )). (c) Total correlations correlations, 
I(Pab)=S (pa) + S (p B ) - S (pab) (see text). Ideal (solid black line) and 
experimental (solid squares with error bars). 



the \d) state to the short-lived P3/2 state (?i/2~7 ns) which, 
due to selection rules, can only spontaneously decay into |0) 
— thereby preparing the mixed state (|aa)(flfl|+|00)(00|)/2. 
Finally, a 729 nm pulse on the |a)<-»|l) transition, followed 
by another on the qubit transition, prepares the target state 
Eq. (1). In summary, the classical correlations inpi are estab- 
lished by first generating quantum correlations, then removing 
all the phase information via an irreversible process on both 
ionic-qubits. 

A local amplitude damping channel, with probability p, is 
realised by moving the corresponding fraction of qubit B's 
population from |1) into the ancillary state \a), by a tightly- 
focused 729 nm laser beam. Next, a circularly polarised 
397 nm laser pulse incoherently transfers the \a) population 
into |0). Ramsey experiments show that coherent superpo- 
sitions in the qubit transition are largely unaffected by this 
397 nm 'optical pumping' pulse. The total experimental time 
is approximately 500 ps, of which 460 ps are required to gen- 
erate the initial state Eq. (1). 

Full state tomography of the generated two-qubit state is 
carried out, for a range of amplitude damping probabilities p. 
The maximum likelihood reconstruction method [21 1 is em- 
ployed, providing the most likely physical state to have pro- 
duced the measurement data. Errors are determined via a 
Monte Carlo simulation of projection noise (due to a finite 
number of measurements) centred around the experimentally 




FIG. 3: Amplitude damping results. Real values of experimen- 
tally reconstructed density matrices for damping probabilities (p) (a) 
0.00 (b) 0.79 and (c) 1.00. Imaginary components are < 0.03. Be- 
low, bar charts show the singular values of the corresponding corre- 
lation matrix. The number of non-zero values gives the rank. The 
fidelities with target states are a) 0.984±0.004 b) 0.992±0.006 and c) 
0.991 ±0.002. The tangle (a measure of two qubit entanglement 1231 ) 
of each state is zero to within error (not shown). 

recorded data. 1000 measurements were taken in each basis 
to convincingly distinguish between zero and non-zero dis- 
cord states. The set of zero discord two-qubit states is of zero 
measure in the total set [22] and therefore any white noise 
(e.g. projection noise) is likely to result in the reconstruction 
of a discordant state. For detailed numerical simulations of 
this effect see Supplemental Material. All quantities are de- 
rived from the reconstructed density matrix, including discord 
via the usual numerical optimisation over projective measure- 
ments 151 IT31 . 

Results from the amplitude damping experiment are sum- 
marised in Fig. [2] Figure [2^ shows the discord observed in 
each direction. The results clearly demonstrate the generation 
of statistically significant amounts of discord Db, while Da re- 
mains zero to within error. In all cases the reconstructed states 
contain no entanglement, quantified by the tangle [23l, to 
within one standard deviation (not shown). Figure|2}5 presents 
the dynamics of the states of qubit B, t ± b, on a cross section 
of the Bloch sphere. The distinguishability is reduced by the 
damping process and consequently, for all 0<p<l, qubit A be- 
comes correlated with non-orthogonal states of qubit B. 

Our results do not imply that the total correlations between 
two systems can be increased via operations on only one sys- 
tem. Fig.pt shows that the total correlations, captured by the 
mutual information /, continually decrease during the damp- 
ing channel. The correct interpretation is that the process con- 
verts some of the preexisting classical correlations, quantified 
by /, into quantum correlations. Indeed, if there are no clas- 
sical correlations present in the initial state, then no discord 
can be generated by operating on only one system ll24ll . Fur- 
thermore, only a very restricted class of discordant states can 
be created by such operations: a set of measure zero in the 
total set Il24ll . That there are fundamentally different kinds of 
discordant two-qubit states, in terms of how they can be gen- 
erated, raises the question of whether there is another way of 
quantifying the correlations in these systems. 

An alternative view on quantum correlations is presented in 
Ref. Il24ll . where in addition to discord, the rank R of the corre- 



lation matrix is considered. Although operations on one sys- 
tem can convert classical to quantum correlations, they cannot 
increase the correlations in terms of R, which is therefore con- 
sidered as an additional quantity of interest. R is calculated in 
the following way. First, the correlation matrix (M) is con- 
structed by writing p in a basis of local Hermitian operators 
(e.g. the Pauli spin operators). R is then the number of non- 
zero singular values of M. Hence, R gives the minimal num- 
ber of bipartite product operators needed to represent a given 
quantum state. Originally, R was introduced as a witness for 
discord 1251 . A state with R-\ is completely uncorrected. If 
R > cLa, where cIa is the dimension of the smallest system, then 
the discord is non-zero. For a system of two qubits d.A-d.B-2. 
The lowest rank of a two-qubit system containing discord is 
2, but in general an R-2 state may or may not contain discord. 
The maximum is R-A, representing a highly correlated state. 
In Ref. [9 1 it is shown that R determines the extent to which a 
state can act as a resource in quantum state transmission. 

Fig. Blpresents a selection of reconstructed states, and their 
singular values, from the amplitude damping experiment. The 
results are largely consistent with an invariant value R=2. 
Small non-zero singular values consistent with a higher rank 
are analysed in the Supplemental Material. Numerical simu- 
lations show that these small but somewhat statistically sig- 
nificant values (larger than zero by less than 2-3 standard de- 
viations) are consistent with the effects of measurement 'pro- 
jection' noise, even with the large numbers of measurements 
employed. The difficulties associated with verifying that sin- 
gular values are strictly zero clearly makes the witness criteria 
for discord very challenging experimentally. 

How then can these strongly correlated high rank dis- 
cordant states be generated? We now show that classi- 
cal noise processes which operate on both qubits are suffi- 
cient to generate states of all ranks. Consider the follow- 
ing: the qubits interact with an environment which causes 
both qubits to suffer identical single-qubit rotations, around 
some axis, by an amount that is not known and fluctu- 
ates from experiment to experiment. Complete dephas- 
ing under this classically correlated noise can be modelled 

by a quantum map ef rf =^; L Kj0)pKl{G) with operators 
Kg(6)=R ft (e)®R ft (0) and Rg(0)=e- m& ' 2 is a single qubit rota- 
tion around a normalized axis vector Ft. The integral over the 
angle 9 generates a dephasing effect between eigenstates of 
the rotation operator. 

This type of noise occurs in any experimental situation with 
fluctuating external fields that couple equally to the qubits. 
In the case of trapped-ion qubits, correlated changes in the 
coupling strength between lasers and ions can occur due to 
intensity or beam-pointing fluctuations, for example. Other 
prime candidates are the fluctuation of ambient electric or 
magnetic fields which remain constant over the closely-spaced 
ion string. 

We investigated the effects of correlated dephasing due to 
ambient magnetic (B) field noise. Although these effects are 
always present, typical experiments require orders of magni- 
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FIG. 4: Correlated dephasing results. Real values of experimentally 
reconstructed density matrices, and corresponding singular values of 
the correlation matrix (underneath), for target states: a) e" d (pi); b) 
p 2 (see text); c) <> 2 ); d) \<fr)=\+,+), where |+>=(|0>+|1»/ V2; e) 
«&(W<0); f) ^ d [K M (7r/2) «* (IM0) K^nliy], for A ■ &=cr y and ft ■ 
&=cr z . The discord D B is a) 0.19±0.03 b) 0.01±0.01 c) 0.188+0.002 
d) 0.011+0.005 e) 0.22+0.03 f) 0.12+0.03. The tangles are all zero 
to within error and the discord D A is equal to D B to within error in 
each case. Imaginary components of density matrices are all < 0.05. 



tude less time than that required for complete dephasing. Here 
we intentionally expose various quantum states for long times. 
Dephasing via B-field fluctuations can be understood as fol- 
lows. The transition frequency of each qubit is determined by 
the local B-field, which is largely identical for both ions due 
to their small separation («5 /mi). B-field noise causes corre- 
lated qubit rotations in the Z basis, i.e. ft ■ cr-cr z , and 6 is pro- 
portional to the B field. Over the many repeated experimental 
runs required to estimate expectation values for quantum state 
reconstruction, the effect of these fluctuations is to smear out 
the reconstructed phases between eigenstates of the cr z opera- 
tors (logical states). (As previously shown [26 1, when applied 
to certain entangled GHZ states this noise process can be ex- 
tremely detrimental, resulting in a fully classical R-2 state. 
Here we show the converse — the process can also generate 
quantum correlations). 

The dephasing channel is applied to the classically corre- 
lated R-2 state p\, by introducing a delay between initialisa- 
tion of the qubits into this state and performing tomographic 
measurements. A delay of 10 ms is found to be sufficient to 
achieve almost complete dephasing, which is over two orders 
of magnitude longer than that required to perform coherent 
laser-driven qubit operations. Experimentally reconstructed 
density matrices and singular values of p\ and e" d (p\) are 
shown in Fig 3. a) and Fig 4. a), respectively for ft-cf — <j z . One 
observes the almost complete loss of the |00)(1 1| coherence el- 
ement of the density matrix |27l . while the |01)(10| coherence 
is largely unaffected since it is in a decoherence-free subspace 
||28l . The measured discord increases from 0.009±0.004 to 
0.19±0.03 and the singular values are consistent with an in- 
crease from an R-2 to an R-3 state. 



In the Supplemental Material we present new theoretical 
results on the conditions under which correlated noise can 
change the correlation rank R of a quantum state and by how 
much. We show that the final rank can be obtained with the aid 
of a simple geometrical picture, expressing the relationship 
between the rotation axis ft of the correlated dephasing and 
two normalized vectors v, w e R 3 which provide all the neces- 
sary information about R-2 states. Specifically, R-2 states 
where the reduced state of each qubit is completely mixed 
(e.g. state p t ) can be written as p = |(I®I + Gfi?-c?(g>vv- <?). 
The final rank depends on the overlap between ft and the vec- 
tors v and w, respectively. Inpi the qubits are correlated in the 
X direction, i.e., v = w — e x and the dephasing rotations are 
conducted around the orthogonal Z direction, ft — e z . In this 
case an R-3 state is generated since ft ■ v — ft ■ w — 0. How- 
ever if ft is neither equal nor orthogonal to v and w then an 
R-4 state can be generated. We demonstrate this by preparing 
the new initial state p2=Kn(7t/S)piKn(n/8y for ft ■ cr-o- y , by 
applying a 729 nm laser pulse Kp,{n/8) to p\. Figure |4jb-c) 
present the results for this initial and final state: the signifi- 
cantly increased 3 rd and 4 th singular values are consistent with 
the conversion from an R-2 to an R-4 state. 

In all cases considered so far classical correlations have 
been present in the initial state. However, we also find that 
completely uncorrected R-\ states (p=pa®Pb) can be con- 
verted to high rank states (see Supplemental Material). The 
conditions for this can be described in terms of the reduced 
Bloch vectors of the two qubits j* 1 and r 8 . If either ft-f 1 ^ or 
ft-r B then R is unchanged by correlated noise. In any other 
case an R-3 state is generated. Figure RJd-e) presents results 
that demonstrate the conversion from an R=l to R=3 state. 
Fig. |4^f) also shows that even conversion from R=l to R-4 
is possible, by effectively employing a sequence of dephasing 
operations. Specifically, as detailed in Fig. HI we allow the the 
R=l state to fully decohere via correlated B-field noise into an 
R=3 state, apply a unitary rotation and then allow the state to 
fully decohere in the same way again. 

Finally we report on the generation of Werner states [29] 
Pw(p)=£# + )(0 + l+(l-.p)I/4, which are a mixture of a maxi- 
mally entangled Bell state </> + =(|00}+|ll»/ V2 and the maxi- 
mally mixed two-qubit state 1/4. These states have rank R-4 
for all p>0 and have been studied in a number of different 
contexts QUI . A range of different Werner states have previ- 
ously been made in light fields [31,32]. For p< 1/3 the Werner 
states are separable, for p> 1/3 they are entangled. In contrast, 
the discord in pw(p), which is symmetric under exchange of 
the particles (Da-Db), runs smoothly from zero at p-0 to a 
maximum at p—\. 

Werner states were made by applying the operation 
MS2- exp(-/ n/4 (TjCTj), where cr I =(cr^+<Ty)/ v2, to the clas- 
sically correlated input state p^.=p|00><00|+(l-p)L The state 
p'w is prepared in the following way: 1) the qubits are en- 
tangled via MS 2(cos- 1 [ yfp\)\00)= VpI 00 >+ V^Pl n > : 2 ) an 
854 nm laser pulse with suitable polarisation causes sponta- 
neous decay of the |11) population, via the P3/2 state, into the 
states \a) and |0) with equal probability — thereby generating 
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FIG. 5: Werner state results, a) Measured discord Da (red circles), 
tangle (blue squares), fidelity with the ideal Werner state (black dia- 
monds). Solid lines show ideal values. Inset shows smooth increase 
in discord, while the tangle (a measure of two qubit entanglement 
1231 ) begins abruptly at p=l/3. The discord D B is equal to D A to 
within error in each case, b) Real values of experimentally recon- 
structed density matrices measured for b) p=0, c) p=0.31 and d) 
p=1.00. Imaginary values are all less than 0.03. Below each ma- 
trix the singular values of the corresponding correlation matrix are 
shown. 



a state equivalent to pS but with the mixed state across the 
|0)<->|a) transition: 3) A 729 nm pulse on the |a)«-»|l) transi- 
tion prepares p'£. Experimental results are presented in Fig.pl 
The observed states achieve a high fidelity with ideal case. 
The discord is found to increase smoothly from zero, while 
the entanglement begins abruptly close to the ideal value of 
p=l/3. 

We have shown that in stark contrast to entanglement, dis- 
cord can be generated between two systems via operations on 
just one system. Not all discordant states can be made this 
way, and the correlation rank provides a way to distinguish 
between different kinds of discordant states. Finally, we have 
shown that noise processes generated by classically fluctuat- 
ing fields are sufficient to generate discordant states with any 
rank, even starting with completely uncorrected states. 
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SUPPLEMENTARY MATERIAL 

Effects of quantum projection noise on tomographic state 
reconstruction 

In this section we present results from numerical simula- 
tions of the effects of measurement noise on experimentally 
reconstructed density matrices. Specifically, we ask the ques- 
tion "if the ideal state is generated in the lab what are the ef- 
fects of taking only a finite set of measurements (i.e. measure- 
ment or projection noise) on the reconstructed density ma- 
trices". The results show that even for very large numbers 
of measurements (a thousand repeated measurements per ob- 
servable) the states reconstructed via maximum likelihood to- 
mography have statistically significant differences to the ideal 
states. 

State tomography requires experimentally estimating the 
probabilities for finding the quantum system(s) to be in var- 
ious states. This must be done by making a finite number of 
repeated measurements. If n trials are made and «,■ instances of 
the outcome i occur then the assigned probability is pj=n.i/n. 
For the expected uncertainty in this value we use the Bino- 
mial proportion confidence interval for one standard deviation 
y/pi(l -pd/n. 

We simulate the effects of a finite number of measurements 
n following a Monte Carlo approach: Taking the probability 
distribution of a perfect state we make a number of copies 
(70) and add in each case noise sampled randomly from bino- 
mial distributions defined by the ideal probabilities and n, the 
number of samples. We then reconstruct the density matrix of 
each noisy data set via maximum likelihood tomography. The 
mean and standard deviation of a range of properties, like the 
discord and entanglement for example, can then be calculated 
and compared to the ideal values. 

As an example we take the ideal fully classically correlated 
and separable state given in equation (1) in the main text: 



Pi 



= ^(|+>A< + I®I+>B< + I + |->A<-| ® |->B<-|) (2) 



This state has zero discord when considering measurements 
on either qubit (Da-Db-0), zero tangle (a measure of two 
qubit entanglement [ 1 1) and correlation matrix singular values 
[CM1,CM2,CM3,CM4]=[1, 1,0,0] giving it a rank R=2. 
Figure [6] presents these properties of the noisy copies of the 
ideal state as a function of the number of simulated measure- 
ments. As one would expect, the more measurements taken 
the closer the properties of the reconstructed states are to those 
of the ideal state. 

Fig.[6h shows that at least 1000 measurements are required 
before the discord and tangle are zero to within one standard 
deviation. In other words, if 800 measurements are taken per 
basis then - even if a perfect separable state was made - one 
would expect to reconstruct a state with non-zero entangle- 
ment beyond one standard deviation 
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FIG. 6: Numerical simulation results showing the effect of taking 
a finite number of measurements on the properties of a tomographic 
reconstruction of pi (Eqn. (1)). See text for more details. The fol- 
lowing properties of the reconstructed density matrices are shown a) 
Discord and Tangle 1 1 1, which are ideally both zero, b) Correlation 
matrix singular values [CM1,CM2,CM3,CM4], which are ideally 
[1, 1, 0, 0] c) Zoom in on CM3 and CM4 in b). Note that discord D A 
is presented and does not differ significantly from D B . 

Fig. |6j)-c shows that the singular values of the correlation 
matrix are also sensitive, yielding a full rank (R-4-) beyond 
one standard deviation even after 1000 measurements. One 
should therefore expect to obtain small but non-zero numbers 
for these values even when generating the ideal state in the lab. 
Note that at some point increasing the measurement number 
(and therefore experimental duration) can lead to other prob- 
lems such as assuring that the identical state is prepared for 
longer times. 

Recall that each point in Fig.|6]is the average outcome of 70 
simulated experiments based on noisy copies of the ideal state. 
FigurelTlshows the full distribution of each of the four correla- 
tion matrix singular values (CM1 - CM4) for 1000 simulated 
measurements of the ideal state. Also shown, as black lines, 
are the experimentally reconstructed maximum likelihood re- 
sults for this state (also taken with 1000 measurements). The 
results show that the experimentally observed non-zero 3 rd 
and 4 th singular values (CM3 & CM4) are entirely consistent 
with those expected for the ideal state. Differences in the 1 st 
and 2 nd singular values (CM1 & CM2) can easily be attributed 
to a slight imbalancing of the, ideally equally weighted, two 
terms in pi (Eqn (1)). This could be caused by small errors in 
the pulse sequence used to generate this state. 



Separable operations can increase the correlation rank 

In this section we show that applying a separable operation 
to a bipartite density matrix can increase the rank of the cor- 
relation matrix [2 1, henceforth called the correlation rank. We 
consider a bipartite Hilbert space 'H - 'Ha ® "J^b and sets of 
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FIG. 7: Numerical simulation results showing distribution of cor- 
relation matrix singular values [CMl,CM2,CM3,CM4] expected 
when performing maximum likelihood tomographic reconstruction 
of the state pi (Eqn. 1) and using 1000 measurements per measure- 
ment basis. The distribution of correlation matrix singular values are 
shown as blue bars. Ideally the correlation matrix singular values are 
[CM 1=1, CM2=1,CM3=0, CM4=0]. Solid black lines are the val- 
ues calculated from the experimentally reconstructed density matrix 
for 1000 measurements. Note that the data shown in this figure was 
used to calculate the last point in Fig. |6b). 



{e B } working on the sub- 



on 



arbitrary quantum operations {i 

spaces 'Ha and 'Hb, respectively. We call an operation e sep 
ti separable if it can be put into the form e sep = Yh=\ p,e^®ef, 
with probabilities p t and P denotes the number of terms in the 
sum. In order to introduce the notation which we will use in 
the following, we first review the definition of the correlation 
matrix. First, we represent the initial state p in terms of arbi- 
trary, fixed bases of Hermitian operators {A,} and {B,} on "Ha 
and fig, respectively: 



4 4 
p = Y J Y j r ij A i ®B j , 

'=1 7=1 



(3) 



where <1a and dg denote the dimensions of 'Ha and 'He- The 
correlation matrix R = (r (J ) can be decomposed into its sin- 
gular value decomposition, using orthogonal matrices U and 
V: R — £/diag(ci, . . . ,c^,0, . . .)V T . Here, c, represent the 
nonzero singular values and L defines the correlation rank. 
Introducing 5, = Yij u ji^j an d Pi - Z; v 'fiBj> where U = (My) 
and V = (vy), we can rewrite p as (2J 



p = ^ j c i S i ®F i . 



(4) 



!=] 



In Ref. |3j it was shown that unilocal operations, defined as 
operations of the form e A ® I, cannot increase the correlation 
rank. First of all, we assure that this also holds for a bilocal 
operation e A ® e B , which requires access to both subsystems. 



p' =(e A <8> e B )p = Z Z Z c i d * e H A k ® B j< 



(5) 



i=\ k=\ 7=1 



with e A (Si) = ZkdikBk and e B (F,) = ZjeijBj- The rank of p' 
is now given by the rank of the matrix F = (/jy), with fy = 
Yif = \ CidikCij which is still limited by the minimum of the rank 
of the three matrices C - diag(ci , . . . , Cl, 0, . . . ), D - (a?,*), 
and E = (gy), hence, by L. Hence, even bilocal operations 
cannot increase the correlation rank. However, if we allow for 
arbitrary separable maps the correlation rank can be increased: 



p' = z pi(4 ® t*)p = z Z Z Z P' c Ak ei A ® b > 

(=1 j=\ k=\ 1=1 



i=I 



(6) 



<l 2 , 



We have introduced ef(Sj) = L^d'-fAk and ef(Fj) = 

d 2 

£ /= fj e'^Bj. Now, the rank of the correlation matrix fy = 
ECi Z/"=i PiCjd' k e' 1 can in principle go up to min{PL, d 2 A ,d 2 B ) 



Correlated dephaslng on two quoits 

We assume a state of two qubits with reduced Bloch vectors 
j* 1 and r B , expressing the reduced density matrices of subsys- 
tems 'Ha and 'Hb as 



Pa 



1 / 
,,= --(1 



^ B -<?). 



(7) 



Here & denotes the vector of Pauli -matrices & - (cr x , cr y , cr z ) T . 
The correlations between the qubits is contained in a real- 
valued matrix B - (fiij), such that the total state can be repre- 
sented in terms of Pauli matrices and the identity in its Fano- 
Form (see, e.g., (4)): 



1 



P = 



: + j] rf o-i 8 I + J^ r? 1 ® <Ti + Z Z B '^ i ® ^ 

(8) 



'=1 ./=! 



The correlation matrix is given by 



R 



1(1 {f) 
4^ B ) 



(9) 



where f denotes the 1x3 matrix r - (r x , r y , r z ). The rank L 
of the correlation matrix R is given by 13 



L = 1 + rk(/3 - f* ® r 8 ) 



(10) 



Here r* ® r 6 denotes the outer product: (j* 1 ® r®)y = r^r^. 
In general, the matrix B can be decomposed into singular val- 
ues, Bij - Yjk=\ v ikdkWjk, with orthogonal matrices V - (vy) 
and W - (w/j). Introducing two sets of orthonormal vectors 



8 



(vk)i - vtk and (vt^).- = Wjk as the columns of V and W, we 
rewrite the state p as 



P = 



: + j* 4 -(?i8)I + I®r B -(?+V c/^v^ • cf ® w k ■ cf 

(11) 



k=\ 



Dephasing by fluctuating rotations 

Single-qubit rotations of a Bloch vector r around an axis 
determined by the normalized vector it are described by the 



operator Ra(8) = e ,9 *°V 2 . Using the relation 

(it ■ cf)(r ■ cf) - (it ■ r)I + i(itx r) ■ cf, (12) 

we obtain 

(it ■ cf)(r ■ cf)(it ■ cf) — (it ■ r)(it ■ c?) — (it X r*X it) • & 

= 2(it-r)(it-cf)-(r-cf). (13) 

In combination with e -"?"°"/ 2 = cos(#/2)I + i s,in(0/2)it ■ cf, this 
allows us to describe an arbitrary rotation of a Bloch vector r 
around it by 



e -i6n-&j2f _ ^ e lW? °V 2 



COs2 \2/ ^"v^) ^'^~ 2i ^ n [^ co& W)^ xi ^' 



+ 2 sin 2 



Q(n-r)it 



(14) 



Averaging uniformly over the angle 8 generates a dephasing 
effect. The new Bloch vector points into the direction it of the 
rotation: 

— f d9e- mcfl2 ?- oV w<?/2 = (it ■ i)ft ■ cf . (15) 
2* Jo 



Correlated dephasing 



A separable map describing correlated dephasing in both 

subsystems is generated by 



1 



In 



e* d (t-cf®rf • cf) = — J d6R H (0)v ■ d-R^O) e R n (0)u- ■ dR'sO) 



2tt Jo 



dee -mm^. foien-m 9 e -m-m$ . & ^n-m 



— -v '• cf ®w ■ cf v ■ cf ®(it ■ w)it ■ cf (it ■ v)it ■ cf® w ■ cf 

1 3 
(vXit)-cf®(wXit)-cf-\ — (it ■ v)it ■ cf® (it ■ w)it ■ cf. 



r 



(16) 



Kraus representation 



with the self-adjoint Kraus operators 



Correlated dephasing can be given in form of a Kraus rep- 
resentation: 



<JP) = \kipK, + l -K 2P K 2 + ^K 3 pK 3 , (17) 



K\ - (-1 ®I + it ■ cf®it ■ cf) 

V2 

K 2 - (I ® I + it ■ cf® it ■ cf) 

V2 

K3 - (-1 ®it-cf + it-cf®I). 

V2 



(18) 



9 



Initial rank-1 states 

We show that L — 1 if and only if p - pa ® ps'- First, if p is 
a product state then obviously L — 1, as p& and pg can be seen 
as elements of operator bases of < K A and 'K B , respectively. 
Conversely, assume that L = 1, then, according to Eq. d5), we 
find that (\,? B ) = citfji) = c 2 (^,^ 2 ) = c 3 (^,/? 3 ) s R 4 
with/? = (p],p2>P3)- From this we conclude that c,- = 1/r 4 
and/?,- = r^r 8 . Inserting this into Eq. <|8p yields 



4 V 



)I+^ 



. + . 



./*• 



cf + 



Z* 4 -, 



./*• 



? ) 



= I(I + ^. (? )®I (I + ^. (? ). 



(19) 



This concludes the proof and implies that rk(/3) < 1 for prod- 
uct states. Moreover, if ik(J3) = 1 then the left- and right- 
singular vectors of /3 are given by the reduced Bloch vectors 
j* 4 and r 8 , respectively and /3 - r^ ® r B , in agreement with 
Eq. ([10]). 

We now determine the rank of the product state after be- 
ing subject to correlated dephasing along ft. The correlation 
matrix of the state after application of the map is given by 



4/R) 



(r 8 • ft)ft 



jr 



4 [ (j* -ft)ft e* d (r* ® r B ) > 



with the rank 



(20) 



rk(«4(£)) = 1 + rk(ef d (^ »?)-(?• rt)(? -ft)it®ft) 

= 1 + rk - (j* 4 ® r 8 - a* 4 ® (r 8 • n)n - (J* 4 ■ ft)ft ® r 8 - (f^ Xft)®(r B Xft) 

+(/*-n)(/* B -n)"®'?)) 
= 1 + rk(i ((t* 4 - (f* 4 • fO/0 ® (? - (i ■ ft)ft) - (f* x ft) is (»" x d))J . 



(21) 



We obtain rk^ ,(R)) - 1 if and only if ft - r^/r* or ft — r B /r B , The matrix ft is given by the outer product [5 - dv®w. 



with r 4 = "Vr 4 " - ^ and r 8 = Vr 8 • r 8 . Otherwise, the final 
rank is 3. 



We will use the following relation: 



( 3 



Initial rank-2 states 

We consider states with maximally mixed reduced density 
matrices, i.e., j* 4 = r B — 0. The rank of the correlation matrix 
is then given by L — ik(J3) + 1. Hence, classical states with 
maximally mixed marginals have ik(J3) < 1 and p 1 - corre- 
sponds to the overall maximally mixed state, which has L = 1. 
Let's assume we have L — 2, hence, only one singular value 
of /3 is nonzero. Such a state has always zero discord as it can 
be written as 



Po = ^fl< 



'. + dv ■ & ® w ■ (?) . 



(22) 



rk 



z 

U=l 



v k ® Wjfe 



rk(VW y ) < min{rkV,rkW}, 



(23) 



where the rank of the matrices V and W is determined by the 
number of linear independent vectors i^ and w/t, respectively. 
If one of the two sets is linearly independent, the rank is given 
by the number of linear independent vectors of the other set. 
Using this it is possible to specify a combination of conditions 
which guarantee the conversion of a classical rank-2 state with 
maximally mixed marginals into a fully correlated rank-4 state 



by correlated dephasing. Correlated dephasing on ( 22 1 yields 
the state 



4m) = 7 



'. -\ — [v • & ® (w — (ft ■ w)ft) ■ & + (ft ■ v)ft ■ & ® (3(ft ■ w)ft — w) ■ & 
-(v x ft) ■ & ® (w x ft) ■ cr]) . 



(24) 



The set {v, (ft ■ v)ft, —v x ft} is linearly independent if and only if ft is neither equal nor orthogonal to v (or put equivalently 



10 



< it ■ v < 1, since both v and it are normalized to one). 
The latter implies that it should be different from all the left- 
singular vectors of /3, even those with singular value zero. Ac- 



cording to Eq. (23 i, the rank of /3 after the correlated dephas- 



ing is now given by the number of linearly independent vec- 
tors in the set {w-(it-w)it, 3(it-w)it-w, wxit}. Again, all three 
vectors are linearly independent if and only if < it ■ w < 1. 
A rank-2 state with maximally mixed marginals is therefore 
converted into a rank-4 state if and only if < it ■ v < 1 and 
< it ■ w < 1 . For states which are symmetric under permuta- 
tion of the two qubits (v = w), this can be expressed directly 
in terms of/3 as the condition 



< it l /3it/d < 1, 



(25) 



where the Hilbert Schmidt norm of yS is given by d - \\/3\\ - 

The above analysis tells us how to create rank-4 states. 
When do we obtain a state of rank 1, 2 or 3? If it — v we 

get 



<4(Po) = t (I ® I - fif(v • w)v- <t® v- <r), (26) 



which has rank 2 if \? ■ w > 0, otherwise rank 1 . The same 
holds for it - w. For it ■ v — the final state 



4M) = 



i 



+ ~ [v • <? ® (w - (n • w)it) ■ & - (y X it) ■ & ® (w X it) ■ &\ 



(27) 



r 



has rank 3 as long as it ■ w < 1 (equivalent to it + w), which is 
also the case for it ■ w — if it ■ v < 1. 



Summary 

For initial rank-1 states p - p A ® p B , with nonzero reduced 
Bloch vectors r^ and r B , correlated dephasing around it yields 
a state of 

L=l if it is equal to either ?** or r 8 , 

L=3 if it differs from both t* 4 and r B . 

For initial rank-2 states with maximally mixed reduced den- 
sity matrices (j* 4 = r B — 0), correlated dephasing around it 
yields a state of 

L=l if it — v (it — w) and it ± w (it ± v), 

L=2 if it — v (it — w) and it jL w (it JL \f), 

L=3 if it + w (it + v) and it ± v (it ± w), 

L=4 if it is neither equal nor orthogonal to and v and w, 

where v and w denote the left- and right-singular vectors of the 
matrix/? = dv® w in Eq. ([8]). These results are summarized in 
Tab. [8] For the special case of v = w, we obtain L — 4 if and 
only if 



<-*£-<!. 

VTr^S 



(28) 



The present analysis can be extended beyond qubit-systems 
based on the definitions of the Fano form, rotations and gen- 
eralized Bloch vectors for higher dimensions [4|. 

it = v? it ■ \$ = < jf ■ v? < 1 



n 



v 



it ■ v = 
< it ■ v < 1 



FIG. 8: Rank of the correlation matrix after application of correlated 
dephasing to a rank-2 state with maximally mixed marginals. 
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